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We study the effect of in-medium hadronic properties in photon nucleus interactions in the context of shadowing as 
well as the dilepton spectrum for incident photon energies in the range 1.1-3 GeV. A reasonable agreement with the 
experimental data for shadowing is obtained in a scenario of downward spectral shift of the hadrons. We show that 
distinguishable features for in-medium changes of the hadronic properties can be observed experimentally through 
the dilepton spectrum by judicious choice of target nuclei and incident energy of photons. 

I. INTRODUCTION 

The study of in-medium properties of hadrons has been a field of great interest for quite some time. The recent 
renewed interest in this area of physics is mainly due to various results available from relativistic heavy ion collision 
experiments performed at the Super Proton Synchrotron at CERN and Relativistic Heavy Ion Collider at Brookhaven 
National Laboratory. Particularly a series of studies by the CERES/NA45, HELIOS-3 and NA50 collaborations at 
SPS with S+Au, S+W, Pb+Au and Pb+Pb collisions have largely been interpreted as the indication of a downward 
shift of the p° mass in the nuclear medium [1,2]. The invariant mass spectrum for dilepton production shows a large 
enhancement at low invariant mass regions. These excess dileptons are thought to originate from the decay of vector 
mesons with reduced mass [3,4] in the medium. 

In a contrasting view, some authors explained this enhancement as a manifestation of the broadening of the p 
spectral function due to its coupling with baryonic resonances [5-7] . It was initially suggested that the drop in p mass 
is mainly governed by the chiral symmetry restoration. But in a theoretical study, combining chiral SU(3) dynamics 
with Vector Meson Dominance (VMD), it was shown that the chiral restoration does not demand a drastic reduction 
of vector meson mass in the nuclear medium [8]. This model also predicted a substantial broadening of p° spectral 
density with only a marginal mass reduction. In ref. [9-11] it was shown that in a non-chiral model like Walecka 
model a drastic reduction in p meson mass may be realized. 

Presently the heavy ion experiments and the corresponding theoretical studies remain inconclusive mainly due to 
the fact that the medium effects here are masked by the complicated dynamics both in the initial as well as the 
final state. Moreover, the huge multiplicities and background makes it difficult to analyze the experimental data 
unambiguously. 

The hadron- nucleus collisions [12], though seem to be less complicated compared to A- A collisions, will also suffer 
from the problem of initial as well as final state interactions. On the other hand, these difficulties are largely overcome 
with the use of photons which do not have the problem of initial state interaction. The high flux photon beam, which 
has compensated for the low interaction cross section with nuclear matter, complimented by wide angle multiparticle 
spectrometers, have made possible a new generation of experiments. 

Detection of dileptons in experiments simplifies the problem further since in contrast to hadronic, dileptonic decay 
modes in nuclear matter are disturbed by the final state interaction. The small cross section for secondary interactions 
of the outgoing electrons with nuclear matter makes the dileptons an ideal probe for studying the reactions inside the 
nucleus. 

The dominant feature of photon interaction above tt production threshold with 0.14 < E 1 (GeV) < 0.5, E 1 being 
the incident photon energy in GeV, is the A resonance production. Above this, in the range 0.5 < E 7 (GeV) < 1 
, photon nucleus interaction can mainly be described by baryon resonance production. The 7 — A cross section per 
nucleon is suppressed compared to the 7 — N cross section in this energy range. Beyond ~ 1 GeV vector meson 
production becomes dominant and VMD is generally used to describe the interaction in this region. 

In the present paper we confine our study in the region 1 < E 7 (GeV) < 3 to understand the role of in medium 
hadrons in shadowing as well as in the dilepton spectrum from 7-A reaction. Our endeavour to study shadowing as 
well as dilepton spectrum is important due to following reasons. 

Firstly, the photo-nuclear data at low photon energies for different nuclei seem to indicate an early onset of shadowing 
[13,14]. This phenomena has been interpreted as a signature for lighter p meson in the nuclear medium [15]. In contrast, 
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the authors in ref. [16] have claimed that the early onset of shadowing can be understood within simple Glauber theory 
[17-20] if one takes the negative real part of the p-N scattering amplitude into account which corresponds to higher 
in-medium p mass. In [21], authors have concluded that the enhancement of shadowing at low energies occur mostly 
due to lighter p meson as well as intermediate n° produced in non-forward scattering. 

But the situation is rather more involved as we have shown in our earlier work [22] . There are several factors which 
might contribute to the mismatch of theory and experiment. Glauber model, which was originally developed for high 
energy scattering, might be improved by taking care of the approximations inherent in it. Furthermore, improved 
estimate of the parameters of the Glauber model as well as multiple scattering, namely cross section/scattering 
amplitudes, might be able to give a better fit to the experimental data. In Ref. [22] we have addressed all these 
related problems to show that a good description of photoabsorption data is not unambiguous and there are some 
serious constraints in the proper understanding of the experimental results. On the other hand, if the dropping vector 
meson mass is the cause of early onset of the shadowing, it should be observed in the dilepton spectrum at low 
energies as well. In fact, the theoretical analysis of the dilepton spectrum is intimately related to the shadowing, as 
the photon- nucleon cross section inside a nucleus has to be multiplied by the effective number of nucleons i.e. A e ff 
(< A, being the mass number) of the nucleus, as obtained from shadowing studies, to obtain the final photon-nucleus 
cross section. Hence a study of the dilepton spectrum might put some additional constraints on the shadowing. The 
dilepton study itself may be less ambiguous as a good description of photon-nucleon data gives a-priori justification 
of the models used. 

Recently, a large reduction of p mass in the nuclear medium has also been reported by TAGX collaboration [23] , 
inferred from the dipion spectrum, p meson is certainly the best candidate to study the medium effects since due to 
its short life time and decay length a large portion of the p° mesons produced will decay inside the nuclear medium. 
But the detection of dipion is disadvantageous as they suffer final state interactions. Though e+e~ branching ratio 
is small, the absence of final state interactions makes the analysis easier. Of course, there will be a large background 
mainly coming from the QED processes like Bcthe-Heitler [24], but these can be eliminated with suitable cuts in the 
final spectrum [25]. 

In the present paper, we have studied, the propagation of the produced vector meson inside the nucleus. The snap 
shots at different times help us to understand the microscopic dynamics in 7-A collisions. It also reveals the fact that 
10 and <p mesons in contrast to the p mostly decay outside the nuclear environment. The path length traversed by the 
vector meson gives a direct correspondence of the dilepton spectrum and the medium density at which it has decayed. 

The energy range considered here is similar to the proposed experiment at CEBAF [26], where incoherent photo- 
production of vector mesons would be studied, with deuterium, carbon, iron and lead as targets. Though, there have 
been theoretical studies of coherent photoproduction earlier, a detailed study with in-medium effects was lacking. The 
present work shows that even in coherent processes, the medium effects would show up in an appropriate kinematic 
window. Furthermore, as the Bethe-Heitler contribution can be suppressed with an energy cut on final dilcptons, the 
results presented here will give a fair idea regarding the quantitative contribution from other coherent processes in 
the presence of in-medium hadrons. 

The paper is organized as follows. In section II, we have described the phenomena of shadowing. Section III is 
devoted to the study of the propagation of vector mesons inside the nucleus and their decay to lepton pairs. Finally 
in section IV we present summary and discussions. 

II. SHADOWING IN PHOTOABSORPTION 

The photonuclear cross section has been found to be suppressed compared to the photon-nucleon interaction for 
photon energies, E y > 1 GeV. This phenomenon is attributed to the phenomena of the nuclear shadowing. Here 
we would like to consider photon energies above 1 GeV where the the dominant process is the production of vector 
mesons in the initial state and the main tool for theoretical description is VMD. 

A. Shadowing in 7-A reaction 

If a photon impinges on a nucleus then, one would naively expect that the nucleus being transparent to the photon 
each nucleon will have equal probability to interact with the beam. In such a case the total cross section for 7-A 
reaction would be given by 

<7 7 A = Aa yN (1) 
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On the other hand, if the photon interacts with nucleus through its hadronic components, then, just like hadron- 
nucleus interaction, the photon's initial interaction will be principally with nucleons in the incoming side of the 
nucleus [27]. The nucleons further along the photon trajectory do not contribute to the total cross section. Hence the 
total cross section is expected to be smaller compared to the previous case, i.e. <j{A) < Ag^-n- This phenomena, 
observed both in hadron-nucleus as well as 7-nucleus interactions, is generally known as shadowing. In case of 7- A 
collision, such a picture would be valid only if the hadronic components of the photon are long lived. This is quantified 
in terms of the coherent length described later. The effect of shadowing can be written as, 

A eff = °lA = 1 + Sa 7 A ^ 

A Aa-yN Acr^N 

where o 1 a = Aa^N+Sc-yA consists of the incoherent scattering of the photon from individual nucleons and a correction 
due to the coherent interaction with several nucleons. In the present section we will discuss the medium effects on 
shadowing in 7-A reactions in the frame work of Glauber formalism [17,18] and multiple scattering [21,28,29] along 
with VMD [30]. 



1. Kinematics and coherence length 

We consider a photon of energy E^r, colliding with a nucleus at rest. The nucleons inside the nucleus move with a 
Fermi momenta pp or energy Ep = \/pp + m*^, which is a function of space co-ordinate through the density n(r). 
The photon energy in the nucleon rest frame then becomes, 

E 1 = lF E lL {l-p F cos6 L ), (3) 

Pp = pp/Ep and 9^ being the angle between incident photon and the Fermi momenta. The total invariant energy s 
can then be written as, 

s = (p 1 +p F ) 2 

= m* N 2 + 2 lF m* N E lL {\~ (i F cos 6 L ), (4) 

where m* N is the effective nucleon mass inside the nucleus. The modification of hadronic masses in nuclear environ- 
ment have been studied in different models [11,31-33]. In the present study we consider two possible approaches 
for the effective mass of the nucleon inside a nucleus: (i) universal scaling scenario (USS) [31] and (ii) Quantum 
Hadrodynamical Model (QHD) [34]. In case of (i) the hadronic masses {run), except those of pseudo-scalar mesons, 
vary with the nuclear density n(r) as 

^ = 1 - 0.2a;, (5) 
m H 

where x = n(r)/n n (r) and n (r) is the normal nuclear matter density (~ 0.15 fm— 3). 

In the QHD model the effective masses of nucleons and vector mesons are calculated using standard techniques of 
thermal field theory [35,10] and parametrized as a function of n(r) as follows: 
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S. = l + T i a j aP. (6) 



For nucleons ai = -0.351277 and a 2 = 0.0766239; in case of p, a x = -1.30966, a 2 = 1.78784, a 3 = -1.17524 and 
a 4 = 0.294456 and finally for u, a x = -0.470454, a 2 = 0.313825 and a 3 = -0.0731274. No medium effects on the <p 
meson is considered as it is expected to be small in QHD [36]. The total invariant energy available for 7./V scattering 
now depends on the position of the participating nucleon through the effective nucleon mass which in turn depends 
on the density, n(r) or the Fermi momentum pp. 

The vector meson produced inside the nucleus will have an effective mass which will depend on the density of the 
nuclear medium as seen by the meson. A direct effect of these mass changes would be reflected in the coherence 
length. The coherence length or the formation length (7 C ) is the length scale of the hadronic component of the photon 
inside the nucleus. In other words, this corresponds to the time scale of the fluctuation between the bare photon 
and the hadronic component of the physical photon. When l c is small, the hadron mediated interaction may become 
indistinguishable from bare photon interaction and there will not be any shadowing. In the present case, l c is a 
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function of the radial distance inside the nucleus. For the vector meson with effective mass my , the coherence length 
is, 



1 



IE. 



£ 7 - Je^ - m* v 2 



m 



2 

*2 ' 



(7) 



where E 7 itself depends on the position of the struck nucleon through Eq. (3) as mentioned before. 

2. Multiple Scattering Formalism 

The formalism for the multiple scattering is based on the optical theorem [28,37]. The nuclear photoabsorption 
cross section in terms of amplitudes (A) of the multiple scattering of the projectile with the nucleons inside the nucleus 
can be written as [29] , 
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where k is the momentum of the photon and n corresponds to the number of nucleons participating in each multiple 
scattering process. The n — l term, i.e., A^ is the amplitude of forward scattering of a photon with one bound 
nucleon and corresponds to the incoherent part in Eq. (2). The nth order scattering amplitude A^ corresponds to 
the process where the incoming photon produces a hadronic state X\ at the first nucleon, X\ propagates freely upto 
second nucleon where a hadronic state X 2 is produced which then propagates freely to the third nucleon. This process 
continues till the hadronic state X n -\ scatters into outgoing photon from nth nucleon. The nucleus is assumed to stay 
in the ground state which means that there is no energy transfer to the n nucleons. Hence the momentum transfer to 
the n t h nucleon is 
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The expression for is, 
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and IIjiq denote the mass and the vacuum self energy of the intermediate hadronic state Xi. M. a b corresponds 
to the invariant amplitude for the process a + n — > b + n and F(qi) is the nuclear form factor, 



Hoi) = jj d\e^- 3 n{x) 



(12) 



where n{x) denotes the nucleon number density. Under eikonal approximation and neglect of the width of the vector 
meson in the intermediate state the above formalism reduces to the Glauber's formula. In the present context, in order 
to understand the role of in-medium hadrons, we modify the Glauber's formula. The photon entering the nucleus at 
an impact parameter b produces a vector meson at position z\. Inside the nucleus, the coherence length l c , in general, 
would be different at z\ (Z c i) and z 2 (^2) as the different densities will yield different masses. The expression for the 
shadowing part of the cross section is then given by [19], 
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where ay = Refvv /I m fvv is the ratio of the real to the imaginary part of the VN forward scattering amplitude [20] . 
<jvn is the V — N scattering cross section [20] and ky is the wave vector of the vector meson. The attenuation of 
the vector meson amplitude is described by the exponential factor. We have included 2-body correlation in the 
two-particle density as [16], (6, z\, z 2 ) — n{b, zi)n(b, z 2 )[l — jo{q c \ z i — z i\)], where q c — 780 MeV and j is the 
spherical Bessel function. The theoretical results of our calculations have been compared with experimental data. 
The extraction of experimental numbers has been discussed in ref. [22]. 

In the present paper we have discussed the results of shadowing for Carbon (C) and lead (Pb) nuclei. Depending 
on the size of the nucleus we have used two different density distributions; for A < 16 the shell model density profile 
of Ref. [38] and for heavier nuclei (^4 > 16) the density profile from ref. [39] has been used. According to eq. (3), E y 
is a function of angle, Ql for non-zero pf- The results which are presented here have been averaged over all possible 
values of 9^. We have found that the effect of Fermi momentum in the kinematics (through eq. (3)) is small. 
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FIG. 1. A e ff/A for C and Pb nuclei as a function of photon energy. We show the results for both multiple scattering 
approach and Glauber model. The dotted, long-dashed and solid lines indicate calculations using Glauber model for vacuum, 
QHD and USS respectively. The circles, dot-dashed and short-dashed lines correspond to the same in the multiple scattering 
approach. 

Fig. (1) shows the variation of A e ff/A with E 1 for C and Pb nuclei. We have used Glauber model and multiple 
scattering approach for the energy ranges considered here. We find that the USS gives a better description of the data 
both in the multiple scattering approach and Glauber model than the scenario with vacuum properties of hadrons. 
In QHD the drop of hadronic masses being larger, in general, the data is underestimated. On an average, our results 
show that the experimental data over the entire range of photon energy under consideration are reasonably well 
reproduced by the downward shift of the spectral function within the USS. 

The results shown here have been found to be sensitive to the pole mass and largely insensitive to the broadening of 
the spectral function. But, more importantly, the results depends crucially on model (Glauber or multiple scattering) 
parameters, and less on the approximations {e.g. eikonal approximation) involved. Since the USS describe the data 
on shadowing reasonably well for the energy range 1 < E~ ( (GeV)< 3, we will evaluate the dileptons within this energy 
range in the USS. For comparison we will give results for both vacuum as well as QHD scenario. 



III. DILEPTONS FROM 7 -A COLLISION 



The photoproduction of vector mesons, specifically p and lu mesons, have been studied before [40,41]. Recently 
dilcpton spectrum have also been studied in photon-nucleus interaction using Boltzmann-Uhcling-Uhlcnbeck (BUU) 
formalism [42]. In this work we consider t channel processes (above resonance region) to study the effects of in- 
medium hadrons on the dilcpton spectrum. The t channel diagrams considered in the present work provides a 
reasonable description of vector meson photoproduction data as we will show below (see also [41]). It is well known 
that to evaluate the lepton pair production from the leptonic decay of the vector meson in a nuclear medium all 
the processes which can either create or annihilate the vector meson under consideration has to be included in the 
imaginary part of its self energy (which is proportional to the in- medium width). 



A. Coupling constants and form factors for t-channel processes 
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FIG. 2. Fcynman diagrams for t-channel processes 



Here we will describe the different processes and the evaluation of the required coupling constants. The reactions 
considered in this paper are shown in Fig. (2). The iso-scalar and iso- vector coupling of electromagnetic current and 
vector fields can be written following the current field identity [43]: 
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jW(I = l) = —JL Pll (14) 
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where m u , and m p are the masses of u>, <fi and p mesons respectively. The couplings g^, g^ and g p are evaluated 
from the e + e~ partial decay widths of the corresponding mesons and found to be 5.03, 13.25 and 17.05 respectively. 

Processes (a)-(d) go through tt and process (e) occurs via er meson exchange. The exchange of vector mesons arc 
not allowed because of the violation of charge conjugation invariance. The r\ exchange is not considered as it will be 
suppressed due to its larger mass and smaller value of coupling compared to pions. For the case of u> in the final 
channel, the neglect of 77 is further justified as ui — > 777 is two orders of magnitude lower than the uj — > 7r°7. Due to 
opposite parity of tt and cr, diagrams (a)-(d), which interfere among themselves, do not interfere with (e) [41]. The 
relative importance of tt and o exchange diagrams can be assessed from the radiative transition of vector mesons. 
The branching ratio of u> — > tt°7 is about an order of magnitude larger than than that of u> — > tt + tt~7 which has an 
upper limit 3.6 xl0~ 3 [44] which implies a dominance of tt exchange for ui in the final channel. In contrast, for p 
meson p — > tt°7 is almost one order of magnitude lower than p — > tt + tt~7 showing the importance of a exchange over 
tt for reactions involving p in the final channel. 

Let us now consider the different interactions needed to evaluate the diagrams. The tt — N and a — N interactions 
are 

CkNN = -igTrNNNj 5 (T.if)N 

CaNN = gaNNNNa (15) 

where g^NN = 13.26 and g„NN — 10.03 [41]. The corresponding vertex form factors are 

„ _ tt(<t) tt(o-) , 

b *NN ~ — T2 -0— (-L6J 
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FIG. 3. Feynman digram for coupling constants evaluation 



The anomalous ujpir interaction [45] is given by 

C Ufm = — ^ d^ v d aP(j tt. (17) 
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FIG. 4. p and ui photo-production cross-sections. The existing data is compared with the theoretical calculations (see [41] 
also). 



The coupling constant g up7T is determined from the decay uj — > 7r° 7 (see Fig. (3a)) assuming a p dominance. The 
decay width is given by 



■ OJ— >7T U 7 



24g2 ml 



,2 \ 3 



mf. 



(18) 



We obtain g W(97r to be 1.92 using r^^ 7r o 7 = 715 keV. The 4>p-K interaction can be written in a similar fashion. By 
using the radiative decay of <fi, g^p^ comes out to be 0.097. 
The cr-exchange interaction for p photoproduction is given by 
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(19) 



The coupling constant g„ pp = 4.6, is determined from p — > ir + n 7 decay width given by (see Fig. (3b)) 
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r P^-7 = 2mp(27r)5 y ^i^ 2 (|mi 2 > (20) 

where E\ and -E2 are the pion energies. The spin averaged matrix element squared is given by 

/i kaP\ - 8?ra dlp P 9l^ 2 (m p -E 1 - E 2 f 

[l l> 3 gl * [2m p (E l+ E 2 )-ml-mir 1 J 

The coupling constant g CT7r7r is obtained as (see Fig. (3c)) 

9mt = ^ (l - 4m»/*£) 1/a = 17.6 (22) 

with m CT = 500 MeV and = 300 MeV. 

The app form factor in monopole form is given by, 

A 2 — m 2 

^VpP A 2 _ „2 ^ 

CPP y 

The results for the differential cross sections of p and uo photoproduction in 7-nucleon collisions along with the 
experimental data for 7 energy 1.4-1.8 GeV are shown in Fig. (4). The cut-off A 7r ( (T ) characterizing the 7r((r)-nucleon 
vertex is 0.7 GeV (1 GeV) [41]. A app = 0.9 GeV is taken from [41]. 



B. Propagation of vector mesons inside nucleus 



The vector mesons produced inside the nucleus take a finite time to travel through the nucleus. The modification 
in their width affects the probability of these mesons to decay inside the nucleus. To have a better understanding of 
this phenomenon, we have studied the propagation of vector mesons with time after they are produced. 

The minimum photon energy required to produce a vector meson of mass M can be calculated as follows. The 
magnitude of the three momentum | | of the vector meson produced in a 7- nucleon collision is, 



P3 



1 



2^ 



(24) 



where A is the triangular function defined in the appendix A. Substituting s — m 2 N + 2m^E 1 in the above equation 
we obtain the threshold for the incident beam energy of the photon as Eq = M 2 /2mjv + M. For the production of p 
meson of mass 770 MeV E = 1.1 GeV. It should be mentioned here that if the mass of the p reduces in the medium 
then this threshold will also reduce. For E n = 1.1 GeV the p mesons will be produced but not u) and (f>. Therefore, any 
interference effects between p and u> will be absent. Moreover, for this incident photon energy the p will be created 
almost at rest so that it decays inside the nucleus. This will lead to a very clear signal for the shift of the p spectral 
function in the medium. 

Now we will describe the methodology used to study vector meson propagation inside the medium. The matrix 
element for 7-N collision inside a nucleus would in general depend on the position r of the nucleon inside the nucleus. 
More specifically, the mass and width of the vector meson are functions of r inside the nucleus and is measured with 
respect to the center of the nucleus. The vector meson produced inside the nucleus at r, would propagate inside the 
nucleus with a velocity v (say) and would decay after traveling a distance given by r' = f + vt in time t. Now, the 
distance traveled by the vector meson would depend on its total decay width in the medium, 



Aot — r + r co n, 



(25) 



where To is the decay width in vacuum and r co ii is the width due to the interaction of the vector meson in the medium. 
The increase in the width of the vector mesons reduce their life time in the medium and hence may allow them to 
decay within the nucleus. The time available to this vector meson to propagate, before it decays into dileptons, would 
be inversely proportional to r to t- So one can now write the distance traveled n = vt as, 
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where the expression of v is, 
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v _\ph}_ / K(r)M 2 

£ 3 "V (s + t-m* 2 (r)) 2 [Z<) 

where all the quantities in Eq. (27) are evaluated at the density of the production point r. The point at which the 
final vector meson decays can now be defined as, 

r' = ^Jr 2 + rf - 2rr t cos 9 (28) 

9 being the angle between the velocity vector v and r. 

The vector meson passing through the nucleus may undergo multiple collisions with the other particles present in 
the medium which leads to a broadening of its width, which can be estimated in following manner [46]. Assuming 
that the vector meson is narrow and medium density is rather low, the intensity of the vector mesons passing through 
the vacuum can be described as, 

m ufe )\ 2 = \m\ 2 e- Tj ^ (29) 

In the nuclear medium this intensity decreases further due to scattering and absorption processes. In the medium Tq 
is replaced by T to t which is given by Eq. (25). So the net reduction would go as exp(— T to ttufe/l)- The increase in 
width due to collisional broadening in the medium is given by, 

Tcoii = n(r)a*wy (30) 

where a* is the total vector meson - nucleon cross section and n(r) is the density of the nuclear medium. The average 
number of collisions the vector meson suffers in the medium can be estimated to be 3 for Pb and 1.6 for C at E y = 1.1 
GeV. For E 1 — 2 GeV these numbers are 3 and 1.6 respectively. These clearly show that number of collisions is not 
very large, even for Pb. In that case, the experimental picture, perhaps will be closer to the no broadening scenario. 
Mass modifications in the medium appears to be the major factor in this case. 



C. Dilepton spectrum 

The formalism discussed in previous sections have been used to evaluate 7- A cross section for different scenarios. 
Here we would like to point out that at low energies some of the baryonic resonances get coupled to the p meson. It 
has been shown in [47] that these couplings make the p spectral function very broad. In the present work this effect 
is taken into account through the width of the spectral function. To avoid double counting one should not separately 
add dilepton production from these resonances as these are already included in the spectral function [48]. 

Now the matrix element for dilepton production (Fig. 2) is recasted so as to incorporate the co-ordinate dependence 
for production and decay points. For a reaction of the type pi + p 2 — > P3 + fci + k 2 this is generically given by, 

lM \2 =A (fci.fc 2 ) [Qpi-fci) 2 + (pi-fc 2 ) 2 ] (P2-P4 + m%(r)) 

1 ' a KW+<W^(r)] [(M2-m2 3 (r0) 2 + m2 3 (r0r2 3 (r')][(t-mL(r))2+rL(r)mL(r)] 1 ) 

where m ex and T ex are the mass and decay width of the exchanged particles ( n and a in the present case). The 
masses of the initial state vector mesons ( m Pl ) produced via VMD and the nucleon mass tun depend on r whereas 
mass m P3 and width r p3 , for the final vector mesons are evaluated at r' . A ex contains all the factors coming from the 
different couplings (see appendix B). The effective masses are evaluated using Eq. (5) for universal scaling and Eq. (6) 
for QHD scenarios. The 4- momenta pi,p 2 ,p 3 , and P4 correspond to the initial photon, initial nucleon, intermediate 
vector meson and final nucleon respectively. The vector meson with momentum p 3 finally decays to lepton pairs of 
momenta k\ and fc 2 via a virtual photon. The differential cross section for such a process can be written as (see 
Appendix for details), 

d(J l+l - 1 l-tlmax 



dM 2 2567^(5,0,™^) 



The detailed forms of the matrix elements M. are given in the appendix along with the expressions as well as allowed 
regions of s, t\ and t 2 . In odcr to incorporate medium effects on the dilepton cross-section from 7 — A collisions 
we convolute the above equation with the density profile of the nucleus concerned. As mentioned earlier, there is 
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shadowing phenomenon in 7 — A reactions which affects the dilepton production. Taking these facts into account we 
obtain dilepton production cross-section per nucleon from 7 — A reactions as 

dM 2 2567T 3 A(s,0,m^) / <2 3 rp(r) 1 j 

The invariant mass distribution of the lepton pairs resulting from the decay of the vector meson inside the nucleus 
is given by 

dM 2 Jp(r)d 3 r { ' 

where P(Ra) is given by, 

P(R A ) = e - M ^r tot /| P 3| ( 35 ) 

is the probability of the vector mesons of mass M to decay outside the nuclear radius Ra- In a similar way one can 
evaluate the dilepton emission rate outside the nucleus by replacing (1 — P{Ra)) by P(Ra) in Eq. (34). 



10" 



> 
c 
3. 



With collisional broadening 
C + y -> C +e + e 

Vacuum 

— uss 

QHD 

Decay outside 



E = 1.1 GeV 

v 



0.3 0.4 0.5 0.6 0.7 0.8 

M (GeV) 

FIG. 5. Invariant mass distribution of lepton pairs from 7-C collisions at E y — 1.1 GeV. The result indicated by vacuum 
corresponds to the mass of the vector meson in vacuum and the in-medium width evaluated by using eq. 25. The curves denoted 
by USS and QHD correspond to the medium dependent masses given by eqs. 5 and 6 respectively. 
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FIG. 6. Same as Fig. 5 for 7-Fe collisions. 
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FIG. 7. Same as Fig. 5 for 7-Pb collisions. 



In Figs. 5, 6,7 we have depicted the the invariant mass distribution of lepton pairs for incident photon energy of 
1.1 GeV for C, Fe and Pb nuclei respectively. da/dM 2 indicates the expected dilepton spectrum if the vector mesons 
decay inside or outside depending on the energy as well as their decay widths. If the meson decays outside, a sharp 
peak is to be observed at the p mass of 770 MeV. In the case of mesons decaying inside, one observes distinguishable 
differences between the vacuum and the in- medium scenarios. At E 1 = 1.1 GeV only the p meson can be produced 
with negligible velocity. In this case the probability for the decay of p inside the nucleus is maximum. Hence it will 
bring the information of the in-medium spectral function of the p most effectively. In the present case the width 
of the p due its interaction with baryonic resonances is given by Eq. (25) and the pole of the spectral function is 
shifted according to Eq. (5) or (6). For all the targets (C, Fe and Pb) considered here, the decay probability of the 
p meson outside the nucleus is much smaller than its decay inside. For incident photon energy close to the threshold 
the heavier the target less is the decay probability outside the nucleus. The effects of the in-medium modification 
of the p is clearly visible through dilepton spectra for the scenarios, i.e. USS and QHD for the entire range of M 
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considered here. 
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FIG. 8. Same as Fig.5 for £ 7 = 1.5 GeV. 
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FIG. 9. Same as Fig.6 for £L = 1.5 GeV. 
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FIG. 10. Same as Fig.7 for £ 7 = 1.5 GeV. 



In Figs. 8, 9 and 10 we show the dilepton emission cross section as a function of the invariant mass, M for E 1 = 1.5 
GeV. At this energy p and u> mesons may be created inside the nucleus with non-zero velocity. E-y = 1.5 GeV is 
below the 4> production threshold. The decay probability of the vector meson outside the nuclear volume increases 
with incident photon energy. For smaller nucleus (carbon, say) the dilepton spectra for the vector mesons decaying 
inside is marginally larger than the spectra originating from the decays outside for M < 0.5 GeV since the vector 
meson with non-zero velocity has a probability to traverse the (small) nucleus before decaying. However, for larger 
nuclei (Fe and Pb) the decay probability inside the nucleus is more. Consequently, the dilepton spectra from the 
vector mesons decaying inside dominates over the one corresponding to the decay outside for invariant mass region 
0.3 < M(GeV)< 0.7 for USS and QHD scenario. 
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FIG. 11. Same as Fig.5 for E„, = 2 GeV. 
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FIG. 12. Same as Fig.7 for E-, = 2 GeV. 

The dilepton spectra for E~ ( — 2 GeV are shown in Figs. 11 and 12 for C and Pb respectively. At this energy all 
the three low mass vector mesons (p, u> and (f>) are produced with finite velocity. The are likely to pass through the 
whole nucleus and decay outside. This is reflected in the invariant mass distribution of dileptons from carbon target, 
where most of the vector mesons decay outside the nuclear volume. In case of Pb target, however, the lepton pairs 
from vector mesons decaying inside dominates in the low invariant mass region. 

In Figs. 13 and 14 we plot the invariant mass distribution of lepton pairs as a function of the path length traveled 
by the vector mesons inside the nucleus for E 1 = 1.1 and 2 GeV respectively. The path length of the vector meson can 

be calculated by using the equation I = -^|J-/r tot . This has been averaged over the Mandelstam variable t (contained 
in ps). I indicates the distance inside the nucleus probed by the vector meson. In view of the fact that at a given 
nuclear density the p is lighter in QHD than USS, it is clear from the results shown in Figs. 13 and 14 that lighter 
the mesons the more it penetrates the target. 
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FIG. 13. The dilepton spectra as a function of average path length inside the Pb nucleus for E y = 1.1 GeV. 
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FIG. 14. The dilepton spectra as a function of average path length inside the Pb nucleus for E 1 = 2 GeV. 

IV. SUMMARY AND DISCUSSIONS 

In this work we have studied the effects of the spectral shift of vector mesons on the photon-nucleus interactions. 
It is found that the experimental data on nuclear shadowing can be well described within the ambit of universal 
scaling scenario proposed by Brown and Rho. The spectral shift of vector mesons in Quantum Hadrodynamical 
model seems to underestimate the data. On the other hand vector mesons with their vacuum properties overestimate 
the data. However, we emphasize that our understanding of the shadowing phenomena vis-a-vis Glauber model may 
be improved through corrections to the approximations inherent in the model as well as by examining the model 
parameters critically. The leading correction to the Glauber model due to deviation from eikonal propagation gives 
rise to a correction ~ {2P^ n R)~ 1 relative to the Glauber scattering amplitude [49]. Here is the centre of mass 
momentum of the nucleus and R is its charge radius. In the present case the correction is rather small for > 1 GeV. 
Moreover, a refinement of the Glauber model (multiple scattering) parameters, e.g., vector meson-nucleon scattering 
amplitude, two nucleon correlation etc. might give a good agreement with the data even with the vacuum properties 
of the hadrons. Hence a better estimate of these quantities is essential for a definitive statement regarding the role of 
medium effects on shadowing in photo-absorption processes. Experimental data with better statistics would certainly 
help us to resolve these uncertainties. 

We have shown that the modification of the hadronic spectral function inside the nucleus drastically change the 
invariant mass distribution of lepton pairs originating from the decays of the vector mesons. To bring out the effect of 
different densities and hence effective in-medium masses and widths, we have explicitly incorporated the co-ordinate 
dependence for production and decay points. The dilepton spectra originating from the vector meson decays show 
strong dependence on the incident photon energy and mass number of the target nuclei. For this purpose we have 
considered incident photon energies of 1.1, 1.5 and 2 GeV on carbon, iron and lead nuclei. It has been demonstrated 
that by tuning the incident photon energy one can create the vector meson inside the nucleus with very small velocity. 
This will facilitate it to decay within the nuclear volume and leptons being electromagnetically interacting particles 
should carry the information of the in-medium properties of the vector mesons very efficiently We have also shown 
that heavier the nucleus more visible are the medium effects. 

It is to be remembered that the "signal" for the modification of the vector mesons properties should be filtered 
out from the background from Dalitz decays and decay of the vector mesons outside the nucleus. By fixing the beam 
energy one can minimize the later contributions. For the estimation of the backgrounds from Dalitz decays one needs 
to know the distributions of the hadrons e.g. n, r], u>, A's etc produced in 7-A collision. The contributions from 
Bethe-Heitler process can be eliminated by suitable cuts in the final spectra. 
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APPENDIX 



A. Differential cross-section for 2^3 particle scattering 

We treat the 2^3 scattering process on the basis of the factorization of the phase space integral into two processes 
2^2 and 1 — > 2. We choose the two-particle intermediate system to be ki+k 2 . The differential cross-section for the 
process pi + p 2 — » p 3 (ki + k 2 ) + Pi is given by [50], 



/' \M~ 2 



J M^E?, (A.1) 



where T is the flux factor given by T = 2 X(s,m 2 ,m 2 ) with X(x,y,z) — x 2 + y 2 + z 2 — 2(xy + yz + zx). The 
Lorentz-invariant phase space factor is given by 

/ dL *» = 1 p$k W (2 *» 4 sHn +»-"-*■-'»)■ < A - 2 > 

In order to write the intermediate state explicitly in the phase space integral, the identity 

1= J ds2 J^ S *(P3 - fc i - k2 ) ( A - 3 ) 
with E 3 = pg + s 2 is used in Eq. (A. 2) to obtain 

/ dLips = j^f J 52 {/ *k iS 5%l + P2 P3 Pi) } 

5 4 {P3 -fa- k 2 ) 



{/ 



d 3 h d 3 k 2 4 



2wi uj 2 

= T^Tg y rfs 2 i? 2 (s, m|, s 2 ) i? 2 (s 2 , , m| 2 ), (A.4) 

where s = (pi +p 2 ) 2 . Let us consider the first two particle phase space R 2 (s,m 2 , s 2 ). Going to CMS: Pi + p 2 = 
one obtains 



-IS 



R 2 (s,mj,s 2 ) = - I 4 EiE3 ${Ei + E 2 - E 4 - E 3 ) 



7T 



2Vi 



J*. 



P4cm / dcosfr 



7 2 4 



2AV2( S 



1 f 

,ml,m 2 2 ) J 



dh, (A.5) 



with ti = [pi - p 3 ) = (P2 - Pi) =m%+m%- 2E 2 E A + 2p 2cm _p icm cos 6 2A . 

The second two- particle phase space integral, in the rest frame of ki + k 2 (quantities denoted by R), can be wriiten 

as 

/d 3 ki 
5(E 3 -ui- lu 2 ) 
4ivibj 2 

= 2, A ' /2fe 8 <-<) | dcosCi (A . 6) 



We define the invariant variable t 2 as follows: 



t 2 = (pi -k 1 f = w 2 + m 2 kl - 2E?u? + 2pf fcf cos0£ pi (A.7) 
so that di 2 = A 1 / 2 (s 2 , m 2 , ii) A 1 ^ 2 (s 2 , m 2 ^ , m 2 2 ) dcos 6j} ipi /2s 2 and we obtain from Eq. (A. 6) 

R 2 (s 2 ,ml,ml) = 2A i/ a( J >ro ; >tl) / ^ ( A - 8 ) 
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Finally, combining Eqs. (A.l), (A. 4), (A. 5) and (A. 8) the dilepton production cross-section is obtained as 

da l+l - \M(M 2 ,t u t 2 )\ 2 



(A.9) 



dM 2 dh dt 2 256tt 3 A(s, 0, m 2 N ) \ 1 / 2 {M 2 , 0, h) 
The limits of t\ and t 2 are given by 

hCt: = 2m% - 1 [(s + m 2 N )(s - M 2 + m 2 N ) T \ 1/2 ( S ,m 2 N ,0)\ 1/2 ( S ,m 2 N ,M 2 )] (A. 10) 

h\™:: = -l(M 2 -t 1 )± [A 1 /2(M 2 ,0,0)A 1 / 2 (M 2 ,i 1 ,0)] (A.ll) 



B. Invariant amplitudes 

The invariant amplitudes for the Feynman diagrams shown in Fig. (2)(a-e) are given by 

M ~ieF^g 1LO g lp g vNN g up7 , e^^pi^e^piXfci + k 2 ) a ui(k 2 )j^vi(k 1 ) u N (p A )^ b u N (p 2 ) . . 

M 2 m^t x -ml\\ml ~ im^T^M 2 - m 2 p + im p T p ] [ ' ' 

M _ -ieF v g lp g^g^ NN g ulm e^ af3 p lp e2(p 1 )(k 1 + k 2 ) a uijk^jpVijk^ u N (p 4 )j 5 u N (p 2 ) 

h M 2 m^[h-ml]\m 2 p -im p Y p }[M 2 -ml+imj: w ] [ ' ' 

M -ieF^g^g^pg^MNg^,^ e^^p^jpiXh + k 2 ) a u;(fc 2 ) 7/3 «;(fc 1 ) u N { PA ) lb u N {p 2 ) 

M^ih-mlttml-im^T^^-mf + impT,] 1 ' > 

M _ -izF lx g lp g 1<i> g 1 ,NNg<t>pT< ^ va ^Pi^l{Pi){ki + k 2 ) a UifahpVijh) u N (p 4 )j 5 u N (p 2 ) 

d ~ M^^-ml^ml-impT^-ml + im^^] [ ' ' 

M _ -^zF a g 2 p g app g aNN e^(pi)[pi • (fci + k 2 )g^ - pf(ki + k^u^k^pv^kx) u N {p4)u N (p 2 ) 
M 2 m^[ti — m 2 + im a T a ][m 2 — im p T p ][M 2 — m 2 + im p T p ] 

The squared invariant amplitude is then 



\M\ = \M a + Mb + M c + Md\ +\M e \ (A.17) 

where 

p2 A 2 2 2 2 

TTT1 — ZJL y-y^yjpy-rrNNyujpir /» ,m 

V la| K + ^r 2 ][(M 2 - m 2 ) 2 + m 2 p T 2 p ] 1 ' * } 

p2 4 2 2 2 2 

' y [m4+m2 r 2 ][(M 2_ m 2 ) 2 +m 2 r 2] \' > 

1 cl k + ^ r 0][( M2 - m w + "OT ' 

[ TO 4+m 2 r2][(M2-m 2 ) 2 +m 2 r 2 ] 1 ' j 
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ReM* a M b = 



p2 A '2 '2 '2 -2 

7r -"-n y^Ljy'ypyTTNNy^piT 



< + mlTl][(M2 - m 2 ) 2 + m^m* + mpT 2 ][(M 2 - + 
(to 2 to 2 + m^mpY^Yp) {(M 2 - to 2 )(M 2 - to 2 ) + m^mpT^Tp} 

- (m 2 m u Y u - mlm p Yp) {(M 2 - ml)m p Y p - (M 2 - m 2 )m w r w } 



(A.22) 



ReM*Mc = 



y^^yi^y^pd^NNd^p-Tg^pTT 



< + mlYl][{M* - m 2 ) 2 + m 2 pY 2 p \ 2 [ m \ + m 2 Y 2 
(to 2 to 2 , + rn^rn^Y^) {(M 2 - to 2 ) 2 + m 2 r 2 } 



(A.23) 



ReM* a M d = 



W y~i^y~i4i9^p9-KN N9upir94>pir 



K + mlr 2 Ji(M2 - m 2 ) 2 + m *r*][m* + m 2 r 2 ][(M 2 - m 2 )^ + m 2 r 2 ] 



(ml 



2_2 



m„ + m. 



.mpY^Yp) { (M 2 - m 2 p )(M 2 - m 2 ,) + mpm^YpY^} 



(m 2 m u T u - mlmpYp) {(M 2 - m 2 )m p Y p - (M 2 - m 2 )m r } 



(A.24) 



^7r ^tt 9joj97<fi9-yp97rNN9ujp7r9(f>pTr 



Hp A//* KA = - - 

b c [m p + m 2 pY 2 p ][(M 2 to 2 ) 2 + m^[m% + "^[(M 2 - to 2 ) 2 + m 2 r 2 ] 

(to 2 ™ 2 , + TO p m r p r ) {{M 2 - to 2 ){M 2 - to 2 ) + m^mpY^Yp) 
(m 2 m < pY < p ~ mlmpYp) { (M 2 - m^wipYp - (M 2 - TO 2 )m w r w } 



(A.25) 



ReM* h M d = 



F-k An 9yuj9f(f>9.yp9 7r l\[l\[gujpTr94>pTr 



H + m 2 pY 2 p ][(M 2 ml) 2 + m 2 r 2 ][(M 2 - to 2 ) 2 + m 2 Y 2 } 

x (M 2 -TO 2 )(M 2 -TO?,) + TO w TO r a) r^ 



(A.26) 



with 



and 



ReM*Md = 



p2 A n 2 n 2 n 2 n 2 
7T ^ 9y(f)9f p9-JT N N 9 ipp-K 



m p + m 2 pYp] \{M 2 + m 2 r 2 ] \m% + m 2 Y 2 } [(M 2 - m 2 ) 2 + m 2 Y 2 } 

{m 2 p m\ + TO p TO r p r ) {(M 2 - to 2 )(M 2 - m|) + TOpTO^TpT^ J 
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\M e \ = 



p2 4 _4 n 2 „2 
o 9 1P 9oNNy<Jpp 



[ mp : + rnpTp][{M*~m^+m 2 pYp] 



A-rr 



2ira \t\ + (t 2 -h + M 2 ) 2 } t x 



M 2 



TO 2 [tl - TO 2 



212 



A„ = 



87ra [t 2 2 + (t 2 -h+ M 2 ) 2 } (4m 2 N - h) 
"M 2 mlliU-mlf+mlYl] 2 



A 2 



to: 



A 2 -h 



m 2 p -m\ 
m p - ti 



F a = 



A 2 - 



A 2 -h 



Kpp - m 2 p 

A 2 p P - h 



(A.27) 
(A.28) 
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